The conductance of interacting nano-wires 
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The conductance of one- dimensional nano-wires of interacting electrons connected to non-interacting 
leads is calculated in the linear response regime. Two different approaches are used: a many-body 
Green function technique and a relation to the persistent current recently proposed based on results 
of the non-interacting case. The conductance is evaluated using the functional renormalization group 
method and the density matrix renormalization group algorithm. Our results give a strong indication 
that the idea of obtaining the conductance from the persistent current holds for interacting systems. 
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The last few years brought fast progress in the experi- 
mental techniques to design and manipulate mesoscopic 
quasi-one-dimensional electron systems. Electron corre- 
lations play a crucial role in one spatial dimension lead- 
ing to Luttinger liquid behavior in the infinite wire limit. 
[1] Furthermore the interplay of electron correlations and 
impurities leads to a drastic effect. Even in the presence 
of a single impurity in the infinite wire limit and at low 
energy scales physical observables behave as if the system 
is split in two chains with open boundary conditions at 
the end points. [2] A deeper theoretical understanding of 
the electron transport through interacting nano-wires is 
thus of great importance in connection with experiments 
and possible applications in nanotechnology. Addition- 
ally it will shed light on the question how electron cor- 
relations and the interplay of correlations and impurities 
influence the physical behavior of finite one-dimensional 
systems. 

In a typical experimental setup the wire - or more gen- 
eral the nano-system - is connected to leads which act as 
reservoirs and the two-terminal conductance is measured. 
The leads are considered to be free of impurities and non- 
interacting. They have different chemical potentials [ii 
and /x r where \ii — /i r is given by the applied voltage. The 
wire might contain impurities and the electrons in the 
wire will later be assumed to be strongly correlated. For 
non-interacting electrons, at temperature T = and in 
linear response the two-terminal conductance is given by 
the Landauer formula G = 2e 2 \T(kp)\ 2 /h with the elec- 
tron charge e and the transmission probability \T{kp)\ 2 
of the wire at the Fermi wave vector kp. [3,4] The factor 
2 is related to the two spin directions. In the following 
we consider spinless fermions and use the dimensionless 
conductance g = Gh/e 2 . The above expression was de- 
rived using scattering theory. For an interacting wire 
at T — and in linear response g can be expressed by 
the one-particle Green function of the wire taken at the 
chemical potential and calculated in the presence of the 
non-interacting semi-infinite leads. [5,6] 

Unfortunately even this Green function cannot be cal- 
culated exactly in most cases of physical interest. Thus 



appropriate approximations are required. We use a func- 
tional renormalization group (RG) method [7,8] which 
we have applied successfully in the context of impurities 
in interacting one-dimensional systems. [9,10] We neglect 
the flow of the two-particle vertex which provides a good 
approximation to the exact result for interactions which 
are not too large (see below). [9,10] This establishes one 
way to approximately determine g. 

Here we exclusively study the zero temperature and 
linear response regime. We focus on the case of an im- 
purity free interacting wire. The generalization to the 
case including impurities is straightforward. The wire 
will be modeled by the lattice model of spinless fermions 
with nearest neighbor interaction U and the leads by a 
spinless tight-binding Hamiltonian. 

Within Luttinger liquid theory the conductance of an 
impurity free interacting wire has been determined using 
an effective field theory and bosonization. [II, 12] Even 
though these works were clearly a step forward in clar- 
ifying the issue that for "perfectly connected" leads the 
conductance is independent of the electron-electron inter- 
action and the length of the wire, i.e. g = 1, in connection 
with experimentally accessible wires the modeling has to 
be doubted. Within the effective field theory it is clear 
what "perfectly connected" means; a sharp jump of the 
Luttinger liquid parameter K = 1 in the leads to K < 1 
in the interacting region. Assuming a spatially varying K 
is an extension of standard Luttinger liquid theory and 
it is not clear how to describe such a situation starting 
from a microscopic model. As later discussed briefly we 
find .9 = 1 for a lattice model in which the interaction 
is turned on infinitely slowly (in space). [13,14] For the 
generic case relevant for experiments g depends on the 
interaction and the length of the wire even if no addi- 
tional one-particle scattering terms are considered at the 
contacts and in the wire. 

Very recently Sushkov suggested to calculate g from 
the persistent current I((f>) of a one-dimensional ring pen- 
etrated by a magnetic flux <fi. [15] The ring of N lattice 
sites is build of two parts: the interacting wire with Nw 
sites, and non-interacting "leads" with Nl sites. The 
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conductance for a fixed Nw is then given by [15,13] 



We consider the Hamiltonian (N = N w + N L ) [20] 
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where Iq{4>) denotes the persistent current of a com- 
pletely non-interacting ring with N — Nw + Nl lattice 
sites. The idea can only be motivated using the concepts 
of the non-interacting case. [16] The leading 1/N behav- 
ior of the persistent current in a non-interacting model 
with N lattice sites and a localized scattering potential 
with transmission probability |T(fc^)| 2 is given by [17] 
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Evaluated at <f> = ir/2 this gives [18] 
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Taking into account the Landauer formula for the non- 
interacting case Eq. (1) follows immediately. If one now 
considers the interacting wire as a localized impurity 
in the otherwise non-interacting ring [13] and assumes 
that the persistent current as well as the conductance 
are given by the same effective transmission probability 
Eq. (1) is also plausible for the interacting case. Using 
non-perturbative methods we give a strong indication, 
that the method indeed works for interacting wires. The 
persistent current can be calculated from the ground- 
state energy Eq(4>) taking the derivative with respect to 
the flux I((f>) — —dEo((j))/d(j). In contrast to the Green 
function method following this suggestion the conduc- 
tance g is accessible using the high precision density ma- 
trix renormalization group (DMRG) algorithm. [19] We 
have earlier calculated the persistent current in interact- 
ing rings using the RG and a complex- valued version of 
the DMRG. [10] This gives a second approach to deter- 
mine g, using either the essentially exact DMRG or the 
approximate RG. 

We show that within the RG both approaches lead to 
the same g and that for not too large interactions the in- 
direct method (via the persistent current) leads to equal 
conductances using RG and DMRG. We thus reach a very 
high consistency using different approaches and methods. 
This gives us confidence that Sushkov's suggestion can 
indeed be used for interacting wires and that the func- 
tional RG approach is applicable also to more complex 
nano-systems. A generalization to wires including im- 
purities is straightforward and within the infinite leads 
approach the conductance in the presence of multichan- 
nel leads can be determined. [14] In contrast to Ref. [13], 
where DMRG is only used for real Hamiltonian matrices 
and the phase sensitivity AE = (L/2)|.Eo(0) — E (tt)\ 
is considered as a measure for the persistent current at 
4> = tt/2 we are able to calculate I (tt/2). 
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in standard second-quantized notation. The hopping ma- 
trix element and the lattice constant are set to one and 
the flux (f> is measured in units of the flux quantum. 
When calculating the persistent current periodic bound- 
ary conditions are used. When using the infinite leads 
method </> is set to zero. The interaction in the wire (lat- 
tice sites 1 to Nw) has to be compensated by an external 
potential. Otherwise the fermions would depopulate the 
interacting region. We here exclusively consider the case 
of a half-filled system. The connection between the leads 
and the wire is considered to be perfect in the sense that 
no local scattering terms, as e.g. a smaller hopping matrix 
element between the sites and 1 (and Nw and Nw + 1) 
or site impurities on the sites and Nw + 1, are included. 
Such terms can be added and the more complex Hamil- 
tonian can still be treated by the methods applied in this 
work. 

In the infinite leads limit and for the model discussed 
here g is given by [5,6] 
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where Gm w ,i is the interacting one-particle Green func- 
tion in the Wannier basis. In Refs. [9] and [10] we de- 
scribe how to determine Green functions using the func- 
tional RG method of a system with interaction on all 
sites. As in these investigations we neglect the flow of 
the two-particle vertex. In the above references it was 
demonstrated that this approximation is valid for inter- 
actions as large as U ss 1. Later it will be shown that 
the same holds here. Crucially within the approximation 
physical observables still show typical Luttinger liquid 
behavior (for U < 2) in the infinite wire limit. [9,10] Thus 
the relevant correlations are included. This has to be con- 
trasted to simple perturbative methods [15,21] which do 
not include Luttinger liquid effects. In the present con- 
text the infinite non-interacting leads can be integrated 
out only leading to additional contributions on the lattice 
sites 1 and Nw- The effective Hamiltonian can then be 
treated using the RG and Gn w ,i(e) can be calculated nu- 
merically. The scattering properties can be understood 
in terms of an effective oscillating one-particle potential. 
This helps to obtain an intuitive understanding of the ob- 
served physics. [9,10] Technical details about this method 
will be presented elsewhere. [14] 
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FIG. 1. Persistent current NI{4>) for N = 64, N w = 12, 
and U = 1. The filled symbols are DMRG data and the open 
symbols RG data. 

In Rcf. [10] we have investigated the persistent current 
of a ring with interaction on all N sites using RG and 
complex-valued DMRG. A generalization to the present 
situation is straightforward. Fig. 1 shows typical curves 
for NI(4>) at even N w obtained from DMRG and RG. 
The parameters are N — 64, Nw — 12, and U = 1. Both 
curves agree quantitatively. This holds for even larger 
interactions (see also Ref. [10]) and for all N and even 
Nw we have investigated. The current for this typical 
set of parameters is far from being sinusoidal. Thus the 
phase sensitivity AEq can obviously not be used directly 
to determine I(tt/2) and g. Only in the limit of a strong 
effective impurity I(<j>) becomes proportional to sin<^>. 
[10] In the present context and for even Nw this limit 
is reached for very large U and Nw- Increasing TV for 
fixed Nw (even) the curves of Fig. 1 only change slightly 
and the infinite leads limit can reliably be extrapolated 
from N w 100. The NI(<f>) extrapolated to the infinite 
lead limit decreases if Nw (even) is increased. This is 
related to the Luttinger liquid scaling of I(<fr). [10,17] For 
odd Nw and all U we have studied the persistent current 
is of almost saw tooth like shape. Increasing N at fixed 
Nw the data seem to converge to the non-interacting 
and impurity free curve NI (<p) — vp(<j) — tt)/tt, with 
< 4> < 2ir and the Fermi velocity vp, [10] which implies 
g = l. To obtain the conductance following the sugges- 
tion of Sushkov we have determined Eq{4>) (and the non- 
interacting groundstate energy Eq (</))) for a few fluxes 
around tt/2 using DMRG and RG. From this I (tt/2) and 
Io(tt/2) were obtained by numerical differentiation. Since 
around <j> = tt/2, Eo(<j>) is a fairly smooth function (see 
also Fig. 1), this procedure leads to reliable results. For 
fixed Nw and U we have done this for N = 16, 32, 64, 
and 128. As an example in Fig. 2 we present DMRG 
data of NI(n/2) for N w = 12 and N w = 13 and dif- 
ferent U as a function of IniV. To extrapolate we have 
fitted the data to a quadratic polynomial in 1/N. For 
even Nw the N — > oo limit of N I (tt/2) can be obtained 
up to high precision. We have also tried different extrap- 



olation schemes to ensure that the asymptotic current 
only very weakly depends on the details of the assumed 
scaling. For odd Nw and U > 3 the finite Nw data still 
change significantly even for N as large as 128 and the 
asymptotic limit obtained by extrapolation is clearly less 
reliable. The DMRG data strongly tend to N I (tt/2) = 1 
which corresponds to g = 1. The same behavior can 
be found for other odd Nw- This suggests g = 1 for all 
odd Nw but a definite statement would require larger N. 
The RG data for NI(n/2) show a quite similar behavior, 
with the difference that using this technique for odd Nw 
the large N data are closer to 1 and clearly extrapolate 
to g = 1. Additionally g was determined within the RG 
using Eq. (5) which already implies infinite leads. 
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FIG. 2. NI{-k/2) for N w = 12 (open symbols), N w = 13 
(filled symbols), different N, and U obtained from DMRG. 

Fig. 3 summarizes our results for Nw — 12 and 
Nw = 13, which are typical also for other Nw- Sev- 
eral conclusions can be drawn from this figure. Within 
the RG the two approaches to determine g lead to results 
which agree quite well. The small deviation between the 
two data sets has mainly two sources, both of them mak- 
ing the persistent current data [Eq. (1)] less precise than 
the infinite leads data [Eq. (5)]: (i) The numerical differ- 
entiation of Eq(4>). (ii) The Nl — > oo extrapolation. The 
estimated total error is smaller than the symbol size. Fig. 
3 shows that Sushkov's suggestion Eq. (1) indeed works 
for an interacting wire. For odd Nw in both approaches 
we find perfect conductance g = 1 for all U. The in- 
finite leads approach does not require an extrapolation 
and solving the flow equations within this approach is 
much faster compared to calculating the persistent cur- 
rent in a TV lattice site system. Thus within the RG and 
also in perturbation theory [21] using Eq. (5) is clearly 
superior to the persistent current approach. The advan- 
tage of the persistent current method is that it enables 
us to determine g using DMRG. For up to U « 1 the RG 
and DMRG data agree quantitatively. For larger U the 
RG data still show the same qualitative behavior as the 
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DMRG results. As discussed in connection with Fig. 2 
the N — * oo extrapolation for odd Nw and large U > 3 
leads to results which are less reliable. The surprising 
even-odd effect, observable also at other Nw, has been 
discussed earlier within the Hubbard model using pertur- 
bation theory. [21] Similar to perturbation theory within 
the RG the effect can be traced back to the combination 
of particle-hole and site inversion symmetry. [14] Using 
the RG the interaction generates a modulation of the 
hopping matrix element between the sites of the inter- 
acting wire. This acts as an effective impurity which for 
odd Nw has a perfect resonance at Uf- For even Nw 
and small U the conductance goes like g = 1 — c(Nw)U 2 
with a c(N\y) > which increases with increasing Nw- 




FIG. 3. Dimensionless conductance g as a function of U for 
N w = 12 and N w = 13. 

For fixed U, g is a decreasing function of Nw (even). In 
the limit of very large Nw we find [14] a power-law sup- 
pression g(U,Nw) ~ N W X y l K ^ known from the prob- 
lem of a single impurity in a Luttingcr liquid. [2] Here 
the "switching on" of the interaction acts as an impurity. 
This behavior cannot be captured by simple perturba- 
tion theory [21] since it requires Luttinger liquid effects 
to be included. In the artificial limit of an interaction 
turned on infinitely slowly (in space) this suppression can 
be avoided leading to g(U,N w ) — 1 also for even N W - 
[13,14] 

In this work, we determine the conductance of an in- 
teracting nano-system (i) by a Green function technique 
and (ii) using the mapping on a persistent current prob- 
lem suggested by Sushkov, [15,13] both at the level of a 
functional RG calculation. The resulting conductances 
are equal and show for interactions U < 1 (iii) excel- 
lent agreement with the conductance obtained follow- 
ing Sushkov's idea and using essentially exact complex- 
valued DMRG calculations of the numerical key quantity, 
the persistent current I(n/2). (i)- (iii) give a very strong 
indication that Sushkov's suggestion, which analytically 
can only be justified for non-interacting nano-systems 



also holds for the interacting case. The conductance can 
be obtained for microscopic models avoiding the mapping 
to an effective field theory using bosonization. As an ap- 
plication we calculate the conductance of an interacting, 
impurity free nano-wire connected to infinite leads. The 
methods used can easily be generalized to (i) more com- 
plex nano-systems with different geometries, (ii) wires 
with interaction and impurities, (iii) more complex con- 
tacts, (iv) multi-channel leads, (v) models with spin, and 
(vi) other filling factors. We thus believe that combining 
the two approaches (infinite lead, persistent current) and 
the two methods (DMRG, RG) gives a very powerful tool 
to investigate transport properties of realistic models for 
complex interacting nano-systems. 
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